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Abstract 

We develop the color dipole gBFKL phenomenology of a diffraction cone 
for photo- and electroproduction 7*A^ — >■ V^A^ of heavy vector mesons 
(charmonium & bottonium) at HERA and in fixed target experiments. 
We predict a substantial shrinkage of the diffraction cone from the 
CERN/FNAL to the HERA range of c.m.s. energy W^. The Q^-controUed 
selectivity to the color dipole size (scanning phenomenon) is shown to 
lead to a decrease of the diffraction slope with (which is supported 
by the available experimental data). We predict an approximate flavor 
independence of the diffraction slope in the scaling variable Q'^ + m'y. For 
diffractive production of the radially excited 25* states (\&', T') the coun- 
terintuitive inequality of diffraction slopes B{2S) < B{1S) is predicted, 
which defies the common wisdom that diffraction slopes are larger for 
reactions with larger size particles. 
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1 Introduction 



Diffractive real and virtual pliotoproduction of ground state ^(15') and radially excited 
V'{2S) vector mesons 

^*P^V{1S)pX{2S)p (1) 

at high c.m.s. energy W = ^/s is an ideal testing ground of ideas on the QCD pomeron 
exchange. The new experimental data on vector meson production coming from the HERA 
and fixed target experiments give a unique insight into how the pomeron exchange evolves 
from the nonperturbative to semiperturbative to perturbative regimes with the increasing 
virtuality of the photon and/or increasing mass my of the produced vector meson and 
have prompted intense theoretical discussions 0, ^, ||, ||, [1^, |ll[ . 

The usual approach to the perturbative QCD (pQCD) pomeron is based on the BFKL 



equation |T^, 13| for the leading-logs (LLs) evolution of the gluon distribution, formulated 
in the scaling approximation of fixed QCD coupling as = const and of infinite gluon correla- 
tion (propagation) radius Rc ( massless gluons); it sums the ladder diagrams with reggeized 
t-channel gluon exchanges. More recently, a novel s-channel approach to the LLs BFKL 
equation has been developed |0 in terms of the color dipole cross section cr(,^, r) (here- 
after r is the color dipole moment, Xg// = (my + Q'^)/{W'^ + Q"^) and ^ = log-^ is the 
rapidity variable). The color dipole approach, to be referred to as the running gBFKL ap- 
proach, is favored because it incorporates consistently the two crucial properties of QCD: 
i) asymptotic freedom (AF), i.e., the running QCD coupling asir) and, ii) the finite prop- 
agation radius Rc of perturbative gluons. AF and the running asir) are an indispensable 
feature of the modern theory of deep inelastic scattering (DIS); without running as{r) it is 
impossible to match the leading-LogQ^ (LLQ^) limit of the gBFKL equation with the con- 
ventional GLDAP equation |]T6| in the overlapping applicability region of the moderately 
small X ([0, |T7|, see also [|r^, [IS], [T^)- The finite Rc is of great importance since the nonper- 
turbative fiuctuations in the QCD vacuum restrict the phase space for the soft perturbative 
(real and virtual) gluons and there is a strong evidence for finite Rc from the lattice QCD 



studies (for the review see [|20[) and hadronic interactions |2^. In the infrared region, 
one also needs to freeze as{r) in order not to run into the Landau pole, asir) < af\ Of 



course, if in our running gBFKL equation [0, |T^] one puts Rc = oo and as = const, then 
the original scaling BFKL equation |T2|, |13[ is recovered [|1^, 

Being formulated in terms of real (rather than reggeized) gluon exchanges, the color 
dipole running gBFKL equation |15| readily incorporates the running as{r). The effect 
of finite Rc can be including by modifying the gluon propagator in the infrared region, for 
instance, introducing the effective gluon mass, fig ~ 1/Rc- Remarkably, in LLs approxima- 
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tion, a finite Rc is consistent with QCD gauge invariance. The freezing of as{r) and the 
gluon correlation radius Rc are the nonperturbative parameters which describe the transition 
from the soft, infrared, to the perturbative, hard, region. The purely perturbative pomeron 
exchange does not exhaust the scattering amplitude and in the practical phenomenology 
of deep inelastic scattering one must add certain soft nonperturbative exchange. It is im- 
portant that the color dipole picture and color dipole factorization for the proton structure 
function and for exclusive diffractive amplitudes do not require the validity of pQCD and 
are viable also for the soft pomeron exchange. The soft pomeron exchange is important only 
for sufficiently large color dipoles, r > R^, and can readily be extracted from the experi- 
mental data on hadronic cross sections and diffractive leptopro duct ion of light vector 
mesons ||^. On the other hand, the pQCD, or hard, pomeron exchange can be related to 
the perturbative gluon structure function of the proton |]14|, p5| . 

Diffractive production of 1^(15*) mesons is particularly interesting because of the so- 
called scanning phenomenon, by which the production amplitude probes the color dipole 
cross section at the dipole size r ~ r^, where 

I ml + g2 



is the scanning radius |, |^ . This scanning property follows from the color dipole factoriza- 
tion for production amplitudes and the shrinkage of the transverse size of the virtual photon 
with and holds beyond the pQCD domain |2^. Varying one can study the transition 
from large, nonperturbative and semiperturbative, dipole size r5 to the perturbative region 
of very short rs <^ -Rc in a very well controlled fashion [T0[. Furthermore, the scanning 
radius rs defines the transverse size of the 7*\^ transition vertex, which contributes to the 
total interaction radius and to the diffraction slope -8(7* V). Motivated by a remarkable 
phenomenological success of such a unified color dipole picture of hard and soft pomeron in 
application to the proton structure function ||2^, ^ and vector meson production p, [1^, in 
this communication we develop the color dipole description of the forward diffraction cone 
5(7* — > V) in exclusive diffractive DIS. We use our early results for the energy dependence 



of the forward cone in color dipole scattering PB[ , obtained from the solution of the running 



gBFKL equation for the diffraction slope Here the crucial point is that breaking of the 
scale invariance by AF, i.e., by running as(r), and finite Rc-, alters dramatically the very 
nature of the BFKL pomeron from a fixed cut in the scaling approximation to a series of 
moving poles for the running gBFKL pomeron (for early quasiclassical analysis see 

also As a result, in the running gBFKL approach one predicts a substantial Regge 

shrinkage of the diffraction slope in the vector meson production, which can be tested at 
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HERA. In this paper we present for the first time the detailed analysis of the depen- 
dence and of the Regge growth of the diffraction slope for the production of charmonium 
and bottonium states. 

Diffractive production of radially excited ^'(25*) mesons will give an additional insight 
into the dipole size dependence of the diffraction slope. Because of the node in the radial 
wave function of the ^'(25") states, there is a strong cancellation between contributions to 
the production amplitude from dipoles r larger than or smaller than the node position r„ 
(the node effect |l|, 0). The resulting strong suppression of diffractive production of 
V'{2S^ vs. V{\S^ has been confirmed experimentally in the J/\E' and \E'' photoproduction 
at HERA and in fixed target experiments Q . It has also interesting manifestations 

in the differential cross sections, which we discuss in the present paper for the first time. 
Because the radius of the V^'(2S') state is larger than the radius of the ground state V(IS'), for 
the diffraction slopes one would naively expect the inequality -8(7* > B{^* —>■ J /'^)- 

However, in this paper we demonstrate that the node effect in conjunction with the color 
dipole factorization predicts a counterintuitive inequality -8(7* '^') < -8(7* J/^)i 
which can be tested at HERA. Because the node effect is sensitive to the form of the dipole 
cross section and its variation with energy, we predict specific energy dependence of the 
V'(2S) /V{1S) production ratios, which also can be tested at HERA. 

The presentation is organized as follows. The subject of the introductory section 2 is 
the color dipole factorization and the determination of the pQCD factorization scales for 
diffractive production. The running gBFKL formalism for the calculation of the color dipole 
scattering matrix and of the diffraction slope is presented in section 3. The decomposition 
of the diffraction slope into the perturbative and nonperturbative components and on what 
physics controls the W"^, flavor and dependence of the diffraction slope is expounded in 
section 4. In Section 5 we discuss in more detail properties of the soft pomeron exchange in 
the color dipole representation. In section 6 we present the salient features of the soft and 
hard exchanges on an example of VN total cross sections. Predictions from the running 
gBFKL dynamics for the forward and t-integrated vector mesons production cross section 
are reported in section 7. We find a good agreement with the low energy data and the data 
from the HERA collider experiments. The subject of sections 8 is predictions for forward 
cone in diffractive production of V{1S) states with special emphasis on flavor symmetry. 
Section 9 is concerned with the node effect in forward production of V{2S) states. The 
summary and some conclusions are presented in section 10. 



4 



2 Introduction into color dipole factorization and pQCD 
factorization scales for diffractive amplitudes 



The Fock state expansion for the hghtcone meson starts with the qq state; the snapshot of the 
relativistic meson as a qq color dipole. The probability amplitude to find the qq color dipole 
r is precisely the qq wave function, \I'(r, z), where z is the fraction of meson's lightcone 
momentum carried by a quark (the Sudakov lightcone variable). The interaction of the 
relativistic color dipole of the dipole moment r with the target nucleon is quantified by the 
energy dependent color dipole cross section, cr(^,r). The effect of higher Fock states qqg... is 
very important at high energy v. To the LLs and/or LL^ approximation it can be reabsorbed 
into the energy (rapidity) dependence of o"(^,r), which is described by the running gBFKL 



equation [0, |T5|. The dipole cross section is flavor independent and provides a unifled 
description of various diffractive processes. 

In the limit of high photon energy i/, the gg-nucleon the scattering matrix Ai becomes 
diagonal in the mixed (r, 2)-representation. This (r, z) diagonalization derives from the 
large longitudinal coherence length, 

Icoh ~ t;^— 2" ' (3) 

and holds if Icoh ^ Rp , where Rp is a size of the target proton. Because the coherence 
length is a purely kinematical scale |Q, the {r,z) diagonalization does not require the 



applicability of pQCD and must hold also for soft pomeron exchange, i.e. even if the dipole 
size r is large. The necessary condition is that the longitudinal scale Isoft for the soft 
pomeron exchange is small, Igoft ^ hoh, which is the case for instance in the dual parton 
string model |^ or different models of exchange by nonperturbative gluons |]5B|, |^, ^ . For 



the phenomenological success of a unifled color dipole picture of vector meson production 
see ilH. 

Taking advantage of the (r*, z) diagonalization of the scattering matrix Ji4, the amplitude 
for real (virtual) photoproduction of vector mesons with the momentum transfer q can be 

represented in the color dipole factorized form 

1 

-M(7* - V,^,Q',q) = {V\M\^*) = jdzj d''f^>\,{r,z)M{i,r,z,q)^,,{r,z). (4) 



Our normalization is such that da/dt\^^Q = |A^p/167r. In Eq. (^), \E'^*(r, z) and \E'y(r, 2;) 
represent the probability amplitudes to flnd the color dipole of size, r, in the photon and 
quarkonium (vector meson), respectively (for the sake of brevity we suppress the spin in- 
dices), and J^{C,,r, z,q) is an amplitude for elastic scattering of the color dipole on the 



target nucleon. The color dipole distribution in (virtual) photons was derived in 1 24, 14 . 
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The color dipole cross section a{^, r) only depends on the dipole size r but not on the 
q-q momentum partition z. Because in the nonrelativistic heavy quarkonium 2; ^ |, at 
small (fin the diffraction cone one can safely neglect the 2;-dependence of Ai and set 2; = |. 
Hereafter we will suppress the argument z. Hereafter either ^ or Xeff resp. x {x is Bjorken 
variable for the inclusive DIS, the straightforward analysis of the relevant Sudakov variables 
gives the relationship Xe// ~ 2x) will be used whenever convenient. 

We focus on calculating the imaginary part of the scattering amplitude for which there 
is a simple representation in terms of the gluon density matrix (see below). The small real 
part can easily be reconstructed from analyticity considerations 



ReM{C,r,q) = '^--^lmM{C,r,q). (5) 

We suppress the discussion of KeAi, which is consistently included in all numerical results. 

The details of calculation of the diffractive amplitude have been presented elsewhere 
1^, |TU|. For the Vqq vertex function we assume the Lorentz structure r'I'7^\E'V^. For the 



s-channel helicity conservation at small q, transverse (T) photons produce the transversely 
polarized vector mesons and the longitudinally polarized (L) photons (to be more precise, 
scalar photons) produce longitudinally polarized vector mesons. One finds 



lmMT{xeff,Q^,q) = ——^ 



q) I jT^^ {mlKoier)<f^{r-z) - [z' + (1 - zy]eK,{er)drcPir, z)} (6) 



2 ^ N,CvV^^^2^ 



where 



lu^Mdxeff.Q .q) (2vr)^ my 

J dha{xeff,r,q) dz Ko{er) \^[ml + z{l - z)my](f){r, z) - d'^cpir, z)^ (7) 

e^ = ml + z{l - z)Q^ , (8) 



«em is the fine structure constant, A^^^ = 3 is the number of colors, Cy = |, |, | 

are the charge-isospin factors for the p°, u^, 0°, J/^P, T production, respectively and Kq^i{x) 
are the modified Bessel functions. The detailed discussion and parameterization of the light- 
cone radial wave function (f){r,z) of the qq Fock state of the vector meson is given in [pi] |. 
For heavy quarkonia one can safely identify the current and constituent quarks. The terms 
cx KQ{er)(f){r, z) and oc eKi{er)dr(p{r, z) for (T), Ko{er)d'^(j){r, z) for (L) correspond to the 
helicity conserving and helicity-fiip transitions in the 7* ^ qq,V qq vertices, respectively. 
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In the nonrelativistic heavy quarkonia, the hehcity flip transitions are the relativistic cor- 
rections, which become important only at large . Eq. (0) corrects a slight mistake in the 
relativistic correction to the amplitude for production of longitudinally polarized photons 
made in 0. The numerical results of Ref. |^ for the J /'^ are only marginally different from 
those to be reported in this paper. 

The representation for a{x,r,q) in terms of the gluon density matrix (see Fig. 1) 

4:71 f cPk 1 ^1^1 

a{x, r,q) = — J — asl^^) [</o(2?^) " Mkr)]J^{x, k + -q, -k + -q) (9) 

where J{){x) is the usual Bessel function. AF dictates that at the gluon-color dipole ver- 
tex, the QCD running coupling must be taken at the largest relevant virtuality, = 
min{fc^, C^r~^}) , where C ~ 1.5 [^] and ensures the numerically similar results of cal- 
culations in both the mixed (r, z) and the momentum representations. The gluon density 
matrix T[x^ k + ^q,—k + is proportional to the imaginary part of the non-forward gluon- 
nucleon scattering amplitude; at g = it equals the unintegrated gluon structure function 
of the nucleon J-'{x, k, —k) = dG{x, k'^)/d\ogk'^ . Eq. generalizes to the non-forward case 
7^ the formula [B^, M for the dipole cross section 



^ , , vrV f , ^,dk^ A\l - jQ{kr)]dG{xeff,k'') 

Because the function /(y) = 4[1 — Jo{y)]/y'^ can qualitatively be approximated by the 
step-function, f{y) ^ 0{Aa- — y) , where A„ ^ |3^, for small r one readily finds 

a{x, r) = -jr^(^s(,r)G{x, ql) , (11) 



where the gluon structure function enters at the pQCD factorization scale ~ ^ fTj, ^ 



39| . For large dipoles, r > -Rc, one can neglect Jo{kr) in the integrand, and the dipole cross 



section saturates, 

/ N 47r^ f ,,o,dk'^ dG(xeff,k'^) 

-(^„„r>R.) = -Jo:s(e)- ^^J^, (12) 

Next, notice that the integrands of (|^),(|^) are smooth at small r and vanish exponentially 
at r > 1/e due to KQ i{er). Because of the behavior cr{x,r) oc in (pT]), the amplitudes 
(iD)® dominated by the contribution from the dipole size r ^ rs given by Eq. (H) - 
the scanning phenomenon p, ^, Q . The scanning property is best quantified in terms of the 
weight functions WT,L{Q'^,r^) defined by 

MrM\€) = 7:-^%^ f '4'-^^^Wr(QW) , (13) 
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(my + Q^)^ my J r"^ 
where in a somewhat abbreviated form (z = T, L, for the exact integrands see Eqs. 

1 

W.{Q\r') = ^r\ml + Q'f j dz^^y^r, z)^ ,,{t, z) , (15) 
^ 

For the 15* mesons to a good approximation the so-defined Wt,l{Q'^ i'^'^) are sharply peaked 
functions of a natural variable y = logr^(Q^ + my). The height and width of the peak in 
^/-distribution do only weakly vary with and the flavor and the peak position defines 



the scanning radius rs ~ At,l/\/Q'^ + my. Consequently, the leading twist terms in the 
expansion over the relevant short- distance parameter r| oc 1/{Q^ + Triy) are of the form 
(here we consider q = 0) 



,„,^,.^^,,oc^_J_G(..„,,i). (17) 

By virtue of (pTj), here the pQCD scale l = tl,t{Q^ + fny) , where the scale parameter 
tt l can be estimated as 



A 



rT.L ^ . (1^ 



T,L 

For the more direct evaluation of the pQCD factorization scales it is convenient to 

substitute ( p!0D into ([T3|) , ([I^) , which then take the form reminiscent of the fc-factorization 
formulas for F2{x, Q^) H H: 

.n. Mri..„.,U- 0, . ,10) 

where 

p, (o^ p] - f "^("^') - -^o^^^)] w (n^ r^) (91) 

OrAQ ' M - y y -^asiQ' + ml.) {krf ^^'^^^ ' " ^ " ^^^^ 

Because of properties of fiy) and of the sharp peaking of Wt,l{Q'^, r'^) Sut r ^ rs, the weight 
functions Qt,l{Q'^, k"^) are similar to the step function, 

eT,L{Q\k^)^e{ql^^-k^), (22) 

and 

.2 



dk^ dG{Xeff,k^)) ^ ,2\\ ril^ 2\ [ 

k"^ d log k"^ 



e.(Q^ e)) = G{x,ff, ql) J ^e.(g^ e)) = G{x,ff, q^)h{Q') , (23) 
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where the factors 

/.xw-) = yy^ „^ — .;„;y y..w-,r-) (24) 

exhibit only a marginal dependence on Q^. 

For small the scale parameters At,l are close to the nonrelativistic estimate A ~ 6, 
which follows from = 3/e with the nonrelativistic choice z = ^. In general, A^^l > 6 
and increase slowly with for heavy quarkonia At,l{'^) ~ 6 at < 100 GeV^ and 

J/^) ~ 6 at = and ~ 7 at = 100 GeV^ which shows the relativistic 

corrections in the charmonium and bottonium electroproduction are small. The corollary of 
the large scanning radius rs and large values of At,l is a very small scale factor tt^l in the 
pQCD factorization scale §: ttA-J/"^) ~ 0.20, rL(p°) ^ 0.15 and tt{p°) ^ 0.07 - 0.10 for 

~ (10 — 100) GeV^, which are substantially smaller than r ^ 0.25 suggested in and 
r ~ 1 suggested in p. Consequently, the moderate values of attainable at HERA do, 
at the best, correspond to the nonperturbative and semiperturbative values of ^^e soft 
contribution to the vector meson production must be substantial and one must be careful 
with the interpretation of the vector meson production data in terms of the gluon structure 
function. The point is that at < ^j/^ the scanning radius is comparable to the radius 
of the J/\l/, has been overlooked in and the formulas of Ref. for the J/\E' production 
amplitudes in terms of the J/"^ wave function at the origin are too crude. Strictly speaking, 
Eqs. (p!9D, pop and (^) were derived for the hard pQCD exchange when rs < Rc and/or 
for for the perturbatively large l- However, because the color dipole factorization is true 
beyond pQCD, one can extend (0) to the soft pomeron and regard this relationship as an 
operational definition of the nonperturbative gluon distribution in the proton. To the same 
extent, Eqs. (pJ]), ( ]20|) and ( ]23|) can serve as a unique basis for extracting the whole gluon 
distribution, perturbative plus nonperturbative, at small x from the experimental data on 
diffractive vector meson electroproduction at HERA. 

The dominance ([T7|) of the longitudinal amplitude at > follows, as a matter 
of fact, from electromagnetic gauge invariance and as such it is true in any reasonable 
model of vector meson production, the familiar vector dominance model (VDM) included. 
The dependence of M.t,l differs drastically from the VDM prediction M-riYDM) oc 
(m^+Q2) 0'tof(pA^), though: instead of ototipN) in (|T6]) one has a{xeff,rs) oc r| oc l/iQ"^ + 
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3 The diffraction cone in the color dipole gBFKL ap- 
proach 



In the familiar impact-parameter representation for amplitude of elastic scattering of the 
color dipole 

lmM{^, r,q) = 2 f £b exp{-iqb)T{^, r, h) , (25) 



(26) 



the diffraction slope B = —2d\oglmA4/dq'^\^^f^ equals 



where 



A(e,r)= I d%b'T{^,f,b) 



(27) 



Then, the generalization of the color dipole factorization formula (^) to the diffraction slope 
of the reaction 7*p — * Vp reads: 

1 

B{l* ^ V,tQ')lmM{^* ^ V,tQ^q = 0) = JdzJ d'rX{tr)%{r, z)^,,{r, z) . (28) 



We sketch here the running gBFKL equation for A(^,r). The running gBFKL 



equation for the energy dependence of the color dipole cross section reads ^ 



/C®a(e,r) 



W{^,Pi) + ^i^,P2)-'y{^,r)] .(29) 



Pi P2 

Here the kernel /C is related to the wave function squared of the color-singlet qqg state with 
the Weizsacker- Williams (WW) soft gluon, in which r is the q-q separation and pi,2 are 
the q-g and q-g separations in the two-dimensional impact parameter plane. The quantity 
S{p) = ^Ggs{p)Ki{pGP)p = ~gs{p)^ pKo{pGP), where K^i^x) is the modified Bessel func- 
tion, describes a Yukawa screened transverse chromoelectric field of the relativistic quark 
and 



Pg 



gs{Ri)Ki{pGpi)— - gs{R2)Ki{fiGp2)— 



\S{p,)-S{p2) 



(30) 



Pi P2 

describes the flux (the modulus of the Poynting vector) of WW gluons in the qqg state. 
The asymptotic freedom of QCD uniquely prescribes the chromoelectric field be computed 
with the running QCD charge gs{r) 



i4:TTas{r) taken at the shortest relevant distance, 
Ri = min{r, pj} in the qqg system. The particular combination of the three color dipole 
cross sections, 



9 

A(t(pi, p2, r) = -[(t(^, Pi) + a(^, pa) - (t(^, r)] , 



(31) 
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which emerges in the r.h.s. of the gBFKL equation, is precisely a change of the color dipole 
cross section for the presence of the WW gluon in the qqg state. 

At short distances, r,pi 2 Rc = ^/pg, the kernel /C does not depend on the infrared 
cutoff Re- The Yukawa cut off of the long range chromoelectric field which has been used 
in Eqs. (^,^) is the simplest phenomenological option. To the LL^ approximation, this 
cutoff is consistent with gauge invariance. If one sacrifices AF putting gs = const and lifts 
the infrared cutoff by letting Rc —>■ oo, one recovers the scale-invariant kernel JC. Both the 
finite Rc and running as break the scale invariance, the detailed discussion of consequences 
is found in |T^, ^7\, ^ ^ . The principal phenomenon is that because of the lack of strong 
logr^ ordering in the BFKL equation there is an intrusion from hard scattering to the regime 
of soft interactions and vice versa, and the effect of the soft region is especially enhanced 
by AF. In the numerical analysis an infrared freezing as{q^) < <y[s^^ = 0.82 has been 
imposed on the three-flavor, one-loop asiq"^) = 47r/[9 log(A;^/A^)] with A = 0.3 GeV. With 
Rc = 0.27 fm, i.e., pc = 0.75 GeV, we found Ajp = 0.4 |T^, the calculation of Regge 
trajectories of subleading pomeron singularities is reported in ||2^, the emerging succesful 
description of the proton strcuture function at small x is published in p6| , p7 |. 

— * 

In the gBFKL equation (^) has been generalized to the profile function r(^,r, 6), 
where the impact parameter b is defined with respect to the center of the dipole: 

^2 



in-" 



2 

d pi Pa 



gs{Ri)K^{pGPi)- - gs{R2)K,{pGP2)- 

Pl P2 



X [r(e, pi, 6 + \p2) + r(e, P2, & + \pi) - r(e, r, h)] . 



(32) 

The calculation of the impact parameter integral (^) reduces Eq. ( ^21) to Eq. (|^). The 
calculation of the moment (|27|) leads to the integral equation for A(^,r). It is convenient 
to separate from the diffraction slope B{^,r) the purely geometrical term y related to the 
elastic form factor of the color dipole of the dipole moment r and to discuss instead of 
A(^, r) the function r]{C,, r) = A(^, r) — |r^o"(.^, r) , which satisfies the inhomogeneous integral 
equation 

2 



gs{Ri)Ki{pGPi)^ - gs{R2)K,{pGP2)- 

Pl P2 

x[v{^,Pi)+ V{^, P2) - r) + ^ (p? + pI- r^) Hi, p^) + a{i, pi)] 

= /C®r7(e,r) + /?(e,r; 



(33) 



where the inhomogeneous term equals 

r) = C® a(e, = ^ / rf'Pi pI 



,Pi 



P2 



gs{Ri)Ki{pGpi) 9siR2)Ki{pGp2) — 

pi P2 
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x(p? + P^-r2)[^(e,P2) + a(e,Pi)] . (34) 

Beacuse the homogeneous piece of equation ( p3D coincides with the gBFKL equation (|29|), 
asymptotically the dipole cross section o"(^,r) and the solution ?7(^,r) of the homogeneous 
Eq. (|3^) have identical energy dependence. Consequently, the solutions of the homogeneous 
Eq. ( PB| ) give the asymptotically constant contribution to the diffraction cone and if cra{^, r) 
is a solution of Eq. (p9D and rja{i,r) is a solution of Eq. (|33D with the diffraction slope 



Ba{C,,r), then rih{C,,r) = ria{^,r) + Ab ■ aa{i,r), where A6 =const, is also a solution of Eq. 
(0) with the diffraction slope Bf){^,r) = Ba{^,r) + Ab . It is the inhomogeneous term, 
/3(^,r), which gives rise to ?7(^,r) oc ^o"(^,r), i.e., to the asymptotic Regge growth of the 
diffraction slope, B{^,r) = B{^Q,r) + 2a^^, and the Regge term 2a^^ does not depend on 
the size of the dipole r. Parametrically, a'jp oc as{Rc)Rl times a small numerical factor. 



With the above specified infrared parameters a'jp ~ 0.072 GeV ^ was found in [28|, for 
slopes of subleading trajectories see |2^ . 



4 The beam, target and exchange decomposition of 
the diffraction slope 

To have more insight into the dipole-size dependence of the diffraction slope, it is useful 
to look at the scattering amplitude (t(^, r, q) in terms of the gluon density matrix. For our 
purposes, it is sufficient to treat the color structure of the proton in terms of the three 
valence (constituent) quarks. Then, as illustrated graphically in Fig. lb, the unintegrated 
density matrix of gluons can be written as 

J^{x, k + \q,-k + ]-q) = - j (fkiT{^, k + Jg, -k + ^q, ki + \q, -ki + ]-q) 

as{k\)\GM^) - G,{h + \q. -h + ig)] , (35) 

where Gi(g^) and 6*2(^1? 1^2) are the single- and two-quark form factors of the proton probed 
by gluons and T(^, k + q, —k + |g, ki + g, —ki + |g) stands for the propagation function of 
two t-channel gluons. In the Born approximation, 

Tii, k + -g, -k + -g, k^ + -g, -ki + -q) = — ^ ^i^- . 36 

^ ^ ^ 2^ \{k + W + ^^l][{k-\q)^ + ^^l] 

Splitting the color dipole vertex function into two pieces, Vd(g, r) = [Jo{\qT) — Jo(A;r)] = 

[Jo(|gr) — 1] + [1 — jQ{kr)], we obtain a useful decomposition 



c^(^, r,q) = — [Jo[^qr) " IJ j «s(/t k + -g, -k + -q) 

+^1 -^as{^')[MkT)-l]Hx,k + -q,-k + -q) (37) 
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Because of the property (|I^) , the second term has the typical logarithmic k"^ integration. It 
comprises the contributions to the q dependence from the target and exchanged gluons. In 
contrast, such a logarithmic /c^ integration is absent in the first term; here the A;^ integration 
converges at finite k"^ ~ Rc^- 
The emerging representation 

-A 47r.^ ,1 , f cfk , 2\n~/ 7 1^ r 1^ 

• y d%r (e, + g, + ^g-, ^1 + ^g, -^1 + ^g)«5(A;D [Gi{q^) - G^ih + ^g, -h + ^g)] , (38) 



nicely illustrates how the three relevant size parameters in the problem give rise to the three 
major components of the diffraction slope. The g dependence coming from the proton vertex 
function Vp{ki, q) = G'i(g^) — 6*2(^1 + |g, —ki + |g) is controlled by the proton size. The g 
dependence coming from the color dipole vertex function Vd = Jo(|gr) — 1 is controlled by the 
color dipole size r. The g dependence coming from T(^, k + ^g, —k + \q, ki + |g, —ki + ^q) 
depends on the effective k'^,kl which contribute to the scattering amplitude and on the 
gluon propagation radius Re- The latter scale remains important even at large k because 
the properties of the running gBFKL pomeron are controlled by interactions at r ~ i?^- In 
the asymptotic BFKL regime, at small x, the ki and k become azimuthally uncorrelated. 

In order to proceed further, one needs a model for G'i(g^) and 6*2 (/ti, ^2)- The radius of 
the proton Rjy probed by the gluon can be different from the charge radius Rch, still Rch 
serves as a useful scale. The two-quark form factor G2{k + |g, —k + |g) is a steep function 
of fc^ and a smoother function of g^ [^. For instance, for the oscillator wave function of 
the 3-quark proton one readily finds 

G2{k + ig-, -k + ig-)] = Gi(^g2)Gi(3P) . (39) 

A straightforward differentiation gives a transparent decomposition of da{^,r,q)/ dq^ into 
the following four terms 



da{^,r,q) 



dcT^'K^,r,q) 

f=0 i=l 



^ rfg2 



16 rd'^k , 2 



g2=0 



3 J k^ ^ 



dq^ 

S^Lr'J d% J d%T{^,k,-k,h,-h)as{kl)[l - G2{h,h)] 



- [1 - Mkr)] J d%as{kl)[l - G2{h,h)] 

+ ^Rl[l - Ukr)] J d%T{i, k, -k, h, -h)as{kl 



g2=0 
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—R%[1 - Mkr)] / d%Ti^, k, -k, h, -h)asikl)G2ih, K 



(40) 



The following properties of T(^, k + |g, —k + ^g, A;i + — A;i + |g) are important in (|40|): 
First, in the infrared regulated QCD it is nonsingular at k"^ = 0, cf. Eq. (^). Second, 
(modulo to the logarithmic scaling violations) its large-fc^ asymptotic is similar to that of 
the Born term (p^), T(^, k + ^q, —k + ki + ^q, —ki + \q) oc . Third, in the Born 

approximation (after the azimuthal averaging) 



dT (e, k + ig, -A? + ig, A^i + ig, -fci + \q) 



dq^ 



52=0 



(41) 



Fourth, the finding of the asymptotic Regge growth of the diffraction slope in p8| implies 
that in the high energy limit ^ — ^ 00 and for all k, ki 

dT (e, k + ig, ~k + ig, h + H -h + |g) 



9g2 



;2=0 



1 



«lp(e - eo) + 0{Ri) T{^,k + -g, -k + -g, k^ + -q, -k^ + -q) . 



(42) 



2^' 2' 2^' 2 
Consider first the decomposition of the diffraction slope for large dipoles, r > R^.. In this 
limit, the amplitude (|38D is dominated by the contribution from k"^ ijlq = R~^ ^ R]^^ so 

Jd^kal{e)J^{^,k,k) , cf. 



that Jo(^'^)) ^1(3/^;^) <^ 1 and can be neglected and cr{C,,r) = 
Eq. (^. Then, the first term in the expansion (^) can be evaluated as 



rfa«(e,r,g-) 



Similarly, 



dq^ 
c?a(3)(e,r,g-) 



rfg2 



g2=0 



■iA6i(T(^,r) = -^rV(^,r) . 



lA63a(e,r) = -^i?Xe,r) 
2 



(43) 



(44) 



The integrand of the fourth term contains the steeply decreasing two-body form factor 
G2{k, —k), which cuts off the integration at k"^ < R%. Consequently, one must distinguish 
between r < i?jv and r > i?jv A simple estimate, which interpolates between these limiting 
cases, is 



dg^ 



2 ^ 24 ' 'r^ + R% 



The bottom line is that A64 ^ A63. Finally, making use of (^2]), the second term in 
can be estimated as 



da(2)(e,r,g-) 



dg2 



^A62cr(e,r) 



«;p(e-eo) + 0(i?c)l ^(e,r). 



(46) 
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At low energy, in the Born approximation, Eq. 
of the resuhing diffraction slope 



If) gives A62 = 2i?^. The salient feature 



cJ ' 



(47) 



is a presence of the geometrical contributions A61 = and A63 = ^R%. 

For large dipoles r >^ Rc one recovers a sort of additive quark model, in which the 
uncorrelated gluonic clouds build up around the beam and target quarks and antiquarks and 
the terms 0{Rl) and 2a'-jp{^ — ^0) describe the familiar Regge growth of diffraction slope for 
the quark-quark scattering. The opposite limit of small dipoles, r ^ R^, is somewhat more 
tricky. In the second and third term, the integration is cut off by [1 — jQ{kr)] and extends 
up to Afj/r"^, precisely as in the dipole cross section ([To|). Consequently, their contributions 



to the derivative (0) are still given by Eq. and Eq. (0), respectively, so that the Regge 
term and the contribution from the target proton size to expansion (|47|) are retained. The 
contribution from the first term, i.e., from the size of the color dipole, changes dramatically 
and will no longer have the geometric form |r^. Indeed, as we discussed following Eq. (|37|), 
the k"^ integration in the first term in (HDl) converges at k"^ "^Rc- Consequently, in this limit 



^2 _ (j'ij^'i Q]2e can factor out ^^(k^) = as{r) from the integrand. This leads to an 



estimate 



^2 ^2 



g2=0 



16 a.,(fi.r<«- 16 



and, after making use of (|TT]), to 



A61 



RlG{i,A^lRl 



(4J 



(49) 



8 G(e,AJr2) 

Similar considerations give an estimate for the contribution to the diffraction slope from the 
fourth term in (^), which is a negligible correction to A61: 

^^^-t2-^-g(i:xa^- 

More comments on A61 are in order. At asymptotically large ^ and/or asymptotically 
small X, the running gBFKL approach predicts the universal x dependence of the gluon 



structure function 17 



G'(a;,Q2) oc 



«5(Q' 



X 



(51) 



where 7 = ^Aip and /5o = 11 — 3^/. Consequently, in the well developed BFKL regime 



A61 will not depend on energy: 



" as (Re) 



(52) 
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However, at moderately small x values, the x dependence of the gluon structure function 
exhibits strong dependence on the factorization scale, the ratio G'(^, 4f )/G'(^, ^) has a 
substantial x dependence and A6i contributes to the energy dependence of the diffraction 
cone. Specifically, it makes the slope of the effective Regge trajecrory a^(e//) substantially 
larger than the true slope of the leading Pomeron trajectory a'jp [^]. 

To summarize, the geometrical contribution to the diffraction slope from the target 
proton size, A63 = |-R|^, persists for all the dipole sizes (the term A64 which is also associated 
with the proton size is negligibly small in all cases). Although the nonperturbative parameter 
i?^ is not calculable from first principles, its contribution to the diffraction slope varies 
neither with energy nor with the dipole size and can eventually be fixed from the accurate 
experimental data. 

5 Soft pomeron and difFractive scattering of large color 
dipoles 

The need for a soft pomeron contribution in addition to the gBFKL dipole cross section 
described previously is brought about by phenomenological considerations. A viable gBFKL 
phenomenology of the rising component of the proton structure function over the whole 
range of studied at HERA (real photoabsorption included) is obtained if one starts with 
the Born dipole cross section 0-3(1^) as a boundary condition for the gBFKL evolution at 
Xq = 0.03 1^ However, such a crB(r) falls short of the interaction strength at r > Re, 



roughly speaking, for the phenomenological value Rc = 0.27 fm one finds ^"^(r > Ifm) ~ 5 
mb, whereas for the description of soft processes one rather needs the dipole cross section 
~ 50 mb at r > 2 fm. Therefore, at r > R^, the above described perturbative gBFKL 
dipole cross section (which hereafter we supply with the subscript "pt") apt{^,r), must 
be complemented by the contribution from the nonperturbative soft pomeron, anpti^,^)- 
Because in all the cases studied the contribution from apt{^, r) exhausts the rise of the total 
cross sections and/or of the proton structure function, in Refs. p6| , |^ we have modeled 
the soft nonperturbative pomeron by the energy independent crnpt{^,r) = CTnptij). For the 
lack of better theoretical and experimental information as well as simplicity, we make the 
simplest possible assumption that the eikonals for the perturbative and soft interactions 
are additive, which to the lowest order amounts to additivity of the dipole cross sections 
= apt{^,r) + anpt{r). 

The direct determination of the total dipole cross section a{^,r), from the experimental 
data on photo- and leptoproduction of vector mesons is reported in |^ and supports the 
flavor independence of cr(.^,r). Other constraints for crnpt{i^) include real photoproduction 
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|T0|, hadronic diffractive scattering ||2^, nuclear shadowing in deep inelastic scattering 
|4^ , diffractive deep inelastic scattering at HERA nuclear attenuation in photo- 



production of light vector mesons and the onset of color transparency in leptoproduction of 
vector mesons and the the proton structure function at moderate and small p7[ . 
All the results are consistent with the form of the dipole cross section suggested in |^ H, 
a convenient parameterization for which is 

-2 1 r _ (t - h.\^ 1 

(53) 



(ynpt\r) = (To 



1 - ^ Aiexp(-'^ 



i=i 



l + ^Aexp(-i^ 



^=1 



with ao = 41.2 mb, = 1.45, A2 = -0.45, ai = 1.30 fm, as = 0.75 fm, = 0.80, 
D2 = 0.36, bi = 0.88 fm, 62 = 2.08 fm, Ci = 0.53 fm, C2 = 1.14 fm. For a somewhat cruder 
fit with A = we find (Tq = 51.6 mb, Ai = 1.82, A2 = -0.82, ai = 1.05 fm, 03 = 0.72 
fm. For small dipoles, r <^ Rc, this cross section is poorly known because it is swamped by 

There isn't anything unusual in the concept of a nonperturbative cross section. The 
conventional gluon structure function of the photon 



Jo 



always contains a contribution from gluons with soft transverse momenta fc^ < Qq ~1 GeV^, 
which persists at all and equals precisely G{x, Qq), the familiar input to the conventional 
GLDAP analysis of the evolution of parton densities. One is perfectly content with the 
strong sensitivity of the GLDAP evolution to this unknown soft input G{x,Ql), which is 
routinely fixed from fits to the experimental data. In the color dipole approach to DIS, our 
soft dipole cross section anpt{r) plays exactly the same role as the gluon (quark) structure 
functions at a soft scale Qq. Furthermore, it is tempting to reinterpret this soft dipole cross 
section cr„pt(r) in terms of the nonperturbative gluon distribution in the spirit of Eq. ([lO|). 
The models of soft scattering via polarization of the nonperturbative QCD vacuum by 
Nachtmann et al. [^, ^ belongs to this category and gives CTnpt{f) very similar to our 

parameterization (0). In the interesting region of r < (1 1.5)fm, the conservative 

estimate of uncertainties in anptij) is 10-20%, the major source of uncertainty being due 
to absorption corrections. For heavy quarkonia the absorption correction are negligible, 
though 0. 

We shall assume the conventional Regge rise of the diffraction slope for the soft pomeron, 
Bnpti^,r) = ABd{r) + AB^ + 2aj^pf.{^ — ^0) , where ABd{r) and ABn stand for the contri- 
bution from the beam dipole and target nucleon size and ^0 = log As a guidance we take 
the experimental data on the pion-nucleon scattering which suggest a^^^ = 0.15 GeV~^ 
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(for the small a^^j descriptions of nucleon-nucleon scattering see ||4§]). A plausible guess for 
the proton size contribution is 

AB^ = A63 = . (54) 

In the energy independent soft exchange for small dipoles ABd{r) is likely to follow the 
geometric law ABd{r) ^ as in Eq. (|43|) . Extension of this law to large dipoles is 
questionable. The large-r saturation of anptij') as parameterized by (|53D is a simplifying 
assumption, what happens at r ^ Ifm is immaterial because even in hadrons the probability 
of finding large dipoles, r ^ Ifm, is negligible. However, the diffraction slope is more 
sensitive to the large dipole contribution. For instance, if scattering of large dipoles of size 
^ ^ Rn is modeled by scattering of thin classical strings off the strongly absorbing target 
nucleon of radius -Rat, then for large dipoles (strings), r > 2i?Ar, one readily finds the profile 
function r(6, r) 9{R — b) + ^9{b — R)9{R+ — b)a,icsmj , which gives the large-r behavior 
o'nptir) ~ 2i?Arr and the tamed rise of the diffraction slope ABdir ^ Ifm) ~ ^'^^ • This 
consideration suggests the interpolation formula of the form 

_ + aRl 

^^^^> - 8 + aR% ' ^^^^ 

where a is a phenomeno logical parameter, a ~ 1. 

Such a taming of the growth of Bd{r) is supported by the phenomenology of ttN scat- 
tering. Take for pion the oscillator wave function and assume that the gluon probed radius 
of the pion equals the charge radius. Then, the contribution from the pion size to the 
diffraction slope for the purely geometrical form (^31) for Bd{r) gives the unacceptably large 

8 (7r|%[<Tpt(xo,r) +or„pt(r)]|7r) 

Taking for the contribution from the proton size ABn the estimate (^) we end up with 
BttN ^ 15 GeV~^, which substantially exceeds the experimental result Bj^Niy = 200 GeV) = 
9.9 ± 0.1 GeV"^ The discrepancy increases further if one adds to the above theoretical 
estimate the Regge term 2a'^p^{^ — ^0) ~ 1 GeV~^ evaluated using the relationship between 
the Xeff and pion energy Xe// ^ £^ . 

What is the origin of this discrepancy? If o"(^,t) was r independent and if the gluon 
probed and charge radii of the pion were identical, then one would find from (|^) the familiar 

AB^{^o) = l{RlH)n^^GeV-\ (56) 



With our parameterization (|53|) , the soft dipole cross section keeps rising at r ~ Ifm and 
for this reason the matrix element (^) is dominated by larger than in the charge radius 
of the pion and we end up with ABt^ larger than an expectation ( ^61) based on the charge 
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radius of the pion. The matrix element can be made smaller and compatible with the 
experiment at the expense a rapid saturation of the soft cross section for large dipoles, 
(^nptif ^ Ifm) ^ iJtoti'^N), when one shall recover the estimate (|56|). This solution must 
be rejected, though, because it would lead to negligible fluctuations of the soft dipole cross 
section in conflict with the experimental data on the diffraction dissociation of pions, which 
require 

{■n\a'^{uQ,r)\Tx) - (7r|(T(z/o, r)|7r)=^ _ ^ ^ ^^^^ 
(7r|a(i/o,r)|7r)2 

An attempt to retain the geometrical law and still agree with the experiment at the expense 
of taking Ai?Ar ~ must be rejected too. We believe that the string model suggested 
taming of Bdij) Eq. ( ^51) is a more acceptable solution. Hereafter we take A_Bjv = = 
4.8GeV"l Then the pion-nucleon diffraction slope is reproduced with reasonable values 
of the parameter a in the formula (|55D: a = 1.2 for a^^^ = O.lGeV^^ and a = 0.9 for 
^'npt — 0.15 GeV~^. Hereafter we shall use the latter set of parameters. The dependence of 
the corresponding diffraction slopes vs. the dipole size r is shown in Fig. 2. For production 
of heavy vector mesons the sensitivity to soft contribution is marginal. 

6 Soft-hard decomposition of total cross sections for 
VN scattering 

We start presentation of our results from an evaluation of the vector meson-nucleon total 
cross section 

^totiVN) = ^ ,2^t-zY I "^'^ + + " ^^]'} ''^''^^f' -^^^^ 

For the parameterization of lightcone wave functions 0(r, z) of vector mesons see |T^. The 



results for Xeff < a;o = 0.03 are shown in Fig. 3. The smaller is the radius of the vector meson 
V the smaller is the total cross section atotiVN), to a crude approximation, atotiVN) oc Ry, 
excepting the radial excitations 

In Fig. 3a we show separately the soft pomeron contribution to atot{VN). For the J/"^ 
the radius is large, -Rj/* ~ 0.4 fm > Rf. = 0.27 fm and the soft contribution is substantial, for 
the T the soft contribution is a small correction to the dominant perturbative contribution. 
At subasymptotic energies, the gBFKL approach predicts steeper rise with energy for smaller 
dipoles, cf. Eqs. (|TT]), the trend which is clearly seen in Fig. 3a. At asymptotic energies the 
contribution from the rising gBFKL cross section takes over for all channels. In it has 



been observed that for the "magic" radius ~ 0.15 fm ~ ^R^ the gBFKL color dipole cross 
section exhibits the precocious asymptotic energy dependence apt{x,r^) oc x~'^^ . Because 
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Rr ~ 0.18 fm is very close to the "magic" radius ta, the predicted energy dependence of 
the perturbative contribution to atotC^N) is very close to oc VT^Ap ^ iyo.8_ 

The case of the \E'' is interesting for its large radius R^,/ ~ 0.8 fm and large soft contribu- 
tion. Because the \E'' and have very close radii, the useful comparison is with atot{4'N). 
For small W atot{^'N) of the present paper is indeed numerically very close to atot{4>^N) 
calculated in 0, but the rise of atoti'^'N) by ~ 50% from W ~ lOGeV to ~ 500 GeV 
is much weaker than the rise of atot{4'N) by almost a factor 2 over the same energy range. 
With our energy independent anpt{r), the rise of atoti'^'N) is entirely due to the perturbative 
gBFKL cross section crpt{C,,r), which rises with energy more steeply at small r. Although 
the and the (p^ have similar mean square radii, because of the node effect the relative 
contribution of small r for the case of "^'N is smaller than for the case of (pN and this 
explains the counterintuitive difference of the energy dependence of the two cross sections. 

7 DifFractive production cross sections for the IS* states 
J/^ and T 

Now we turn to the vector meson production. The strong point about color dipole factoriza- 
tion equations (|ll), (|l4|) , (^), ( ^0|) is that apart from the trivial factors Cy and Cy jmy the 



production amplitudes are flavor independent when considered as a function of the scanning 
radius rs and/or + i, |, |, ||, . To this end, Eqs. (IT1),(0), dH), © represent 



the leading twist terms and the correct twist expansion goes in powers of + m^) 

rather than in powers of For instance, in |]T0| we have shown how the ratio of the J/\E' 



and p production cross sections becomes remarkably constant when the two cross sections 
are taken at equal + irty in contrast to a variation by about three orders in magnitude 
when the two cross sections are compared at equal Q^. For this reason we strongly advocate 
the presentation of the experimental data as a function of the flavor-symmetry restoring 
variable + my rather than and whenever appropriate we present our results in terms 
of this scaling variable. 

The soft/hard decomposition of production amplitudes depends on the relationship be- 
tween Ts and Rc- The hard contribution dominates at <Rc, i.e., if 

Q^ + ml>^r^30GeV\ (59) 

which holds better for the heavier vector mesons and the larger Q^. Our phenomenological 
soft interaction, as well as other models for the soft pomeron [^, |ll], extends well into 
1^ ^ Rc- Arguably, with better understanding of the perturbative gBFKL amplitude, one 
can eventually use the vector meson production for better fixing the effect of soft interactions 
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at short distances. In Fig. 4 we show our decomposition of the production amphtudes into 
the hard (perturbative) and soft contribution as a function of for different vector mesons 
at the typical HERA energy W = 150 GeV. Because for the hghter mesons the pQCD scale 
parameter is smaller, tl{p^) < tl{J/^) < ti(T), the soft contribution is somewhat larger 
for the lighter quarkonia. 

Regarding a comparison with the experimental data, the most straightforward theoreti- 
cal predictions are for the forward production and we calculate da/dt\t=o and B{t = 0). The 
experimental determination of these quantities requires extrapolations of da/dt to t = 0, 
which is not always possible and one often reports the t-integrated production cross sections. 
The principal lesson from the high precision vr^A^ scattering experiments is that the diffrac- 
tion slope B{t) depends strongly on the region of t and for average (t) ~ 0.1-0.2 GeV^ which 
dominate the integrated total cross section, the diffraction slope is smaller than at t = by 
~ 1 GeV~^ ESi. We take these irN scattering data for the guidance, and for more direct 



comparison with the presently available experimental data instead of the directly calculated 
B{t = 0) in all the cases we report 

B = B{t = 0) -IGeV""^ (60) 

which we also use for the evaluation of the t-integrated production cross section from the 
theoretically calculated da/dt\t=o- 

1 c/(t(7* ^ V) 



a(7* V) 



B dt 



(61) 
t=o 



The uncertainties in the value of B and with the evaluation ( pT| ) presumably do not exceed 
10% and can be reduced when more accurate data will become available. 

We start presentation of our results and comparison with the available experimental data 
with real photoproduction of the J/\l/ in Fig. 5. The agreement with experimental data from 



the fixed target experiments (EMC |9|; E516 E401 0; E687 @) and from the HERA 



experiments (ZEUS HI |^ |^, ^) is good in the both magnitude and energy 

dependence of the cross section. For the J/\l/ and = the scanning radius is still large, 
rs ~ 0.4 fm and Fig. 4 shows that at smaller energies W ~ 15 GeV the soft contribution 
makes ~ 50% of the photoproduction amplitude. Still, it is smaller than in atotiJ/"^ N) and 
(7(7 — s> J/^E') rises much faster that af^f.lJ/'^ N), which is one of examples of the failure of 
the vector dominance model for processes with heavy quarkonia. Recall that VDM predicts 
a(7- J/^)ocaL(J/^Ar). 

The relationship (|T^ (and also (|^)) is to a large extent the model independent one 
and predicts the dominance of ai at large Q^. It is convenient to present the results for 
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i? = o"i/o"r in the form of the ratio 



R 



LT 



(62) 

t=o 



Q2 darir-^V) 

shown in Fig. 6. The point made in ^ and in a somewhat different form repeated in 
IP is that compared to A^^ the transverse amphtude Aix receives larger contribution from 
large r asymmetric end-point configurations with z{l — -C 1. For this reason Rlt < 1 
and is steadily decreasing with Q^. The steeper rise of (Tpt{x, r) at smaller r makes the end- 
point contributions less important at higher energies and Rlt rises with energy, although 
very weakly. The above predictions for R = dai/daT must be tested at t = 0, the present 
experimental data on R are for the t integrated cross sections. In it has been argued that 
at large t rather o"t ^ o'l, so that the ratio R measured experimentally for the t integrated 
cross sections can be somewhat smaller than our predictions for t = 0. The calculation of 
the t-dependence of Rlt is an interesting subject on its own and goes beyond the scope of 
the present analysis. 

In the typical HERA kinematics the longitudinal polarization of the virtual photon e ~ 1 
and as our predictions for polarization unseparated production cross section we present 
(j(7* ^ V) = o-xiY ^ V) + a^i^Y ^ V)- In Fig- 7 we present our predictions for the 
J/\E' and T production. The short hand representation of the same results is in terms of 
the exponent of the energy dependence of the t- integrated a (7* ^ V) oc = W*^^f^ 
and/or da/dt\t=Q oc = W'^^'^ff. The exponent 6 for the t-integrated cross section is 
slightly smaller because of the shrinkage of the diffraction cone. The effective intercept 
Agjj depends on the range of W the fit is made (the more detailed discussion of this issue 
is found in 0), in Fig. 8 we present our evaluations for W = 100 GeV. For the sake of 
completeness, we show on the same plot A^ff for light vector mesons evaluated from cross 
sections calculated in Ref. [|10]. Slight departures from exact flavor symmetry are due to 
slight differences in the pQCD scale factors t{V) for different vector mesons. The predicted 
downwards turn of Ag// at very large is due to the increase of Xg// at fixed W. The 
real photoproduction of T offers one of the best determinations of the intercept Ajp of the 
gBFKL pomeron, because in this case one has the the magic scanning radius rs ~ |-Rc and 
we indeed find ^eff ~ Ajp = 0.4. The usual fits to the experimental data are of the form 
a{^*p Vp) oc = W^^^ff. The evaluated value of 5 ~ 0.9 from Fig. 5 in the range 
AO <W < lAOGeV is in good agreement with the value 6 = 0.92 ± 0.14(stat.) ±O.W{syst.) 
extracted from the data on elastic J/\I' photoproduction Analogous estimation of 

6 ~ 0.82 from Fig. 5 in the range 30 < ly < 2A0GeV is in good agreement with the value 



6 = 0.77 ± 0.13 presented in The recent HI data on elastic virtual photoproduction 
of J/^ H reported the values of 5 = 1.07 ± 0.57 at Q"^ = 3.7 GeV^ and 1.22 ± 0.52 at 
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= 16 GeV^ in the energy range 40 < W < 150 GeV, which correspond with our results 
6 = 0.98 and 6 = 1.15, respectively. 

The salient features of the dependence are best seen when cross sections are plotted 
as a function of the flavor symmetry restoring variable + my, when the J/\E' and T 
production cross section exhibit very similar dependence (Fig. 9). With Rlt ~ 1 the theory 
predicts 

If one fits ([63| ) to the (Q^+m^)"" behavior and neglects the dependence coming from the 
gluons structure function, then n ~ 3. The effect of the gluon structure function is twofold. 
At a fixed Xg//, i.e., when energy changes with according to W"^ = {Q"^ + my)/xeff, the 
gluon structure function rises with Q^, which lowers the fitted exponent n: n <3. On the 
other hand, experimentally one usually studies the dependence at fixed energy W, when 
Xeff = {Q^ + my)/W'^ increases with Q"^. Because the gluon structure decreases towards 
large x, this induced dependence enhances the exponent n. The exponent n depends on 
the range of the fit is performed in. For instance, in the J/\l' production at a typical 
HERA energy W = 100 GeV we predict n ^ 2.8 for the semiperturbative region of < 10 
GeV^ and n ^ 3.2 if the fit is performed in the region of 15 < < 100 GeV^. We recall that 
for the p° production we found n ^ 2.4 < 10 GeV^ and n ^ 3.2 for 15 < < loo GeV^ 
|T0| . The results for the T are similar to the large-Q^ result for the J/\I'. The departures 
from the exact fiavor symmetry due to Rlt 7^ 1 and slight fiavor dependence of the pQCD 
scale t{V) are marginal for all the practical purposes. 

The experimental data on virtual photoproduction of charmonium states are still scanty 
and there are no data yet on the photoproduction of bottonium. In Fig. 10 we present 
a summary of the experimental data on the J/\E' production from the fixed target EMC 
experiment [4^ and the ZEUS 59| and HI |]56| , |58| , |60[] experiments at HERA. The theo- 
retical results are for ly = 15 GeV appropriate for the EMC experiment (the dashed curve) 
and for W = 100 GeV appropriate for the HERA experiments, there is a reasonable agree- 
ment between the theory and experiment. One of the outstanding experimental problems 
at large is a separation of elastic reaction 7*p —>■ V + p from the inelastic background 
7*p ^ V + X and the low energy EMC data are well known to have been plagued by 
the inelastic background. The contribution from inelastic background to the experimental 
cross section may be the reason why we underestimate the experimental data. One more 
argument on favor of this point will be presented in the discussion below of the diffraction 
slope. 
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8 Diffraction cone for the V{1S) states 



Evidently, the contribution to the diffraction slope from the 'y*-V transition vertex decreases 
with the decreasing scanning radius rs, i.e., with rising Q"^ At fixed energy W the value 
of Xeff rises and the rapidity ^ decreases which also diminishes the diffraction slope because 
the Regge term becomes smaller which is an important component of the dependence at 
fixed W. In this section we report evaluations of the diffraction slope based on Eq. (pS]). We 
use the results of Ref. p8[ for the energy and dipole size dependence of B{^,r) for gBFKL 
color dipole amplitude. For the soft pomeron contribution, we use the parameterizations 
(^) and (|5BD. According to Fig. 4, the nonperturbative contribution to the J/\E' and T 
production amplitudes is small and our results for the diffraction slope are insensitive to 
the soft pomeron effects. Our definition of the diffraction slope is Eq. (^) in section 7 and 
is meant to correspond to the experimentally measured slope B(t) at t ~ (0.1-0.15) GeV^^. 

In Fig. 11 we present the predicted energy dependence of the diffraction slope for real 
photoproduction of the J/\I' and T. As it was shown in [^, at subasymptotic energies 
the diffraction slope for the gBFKL amplitude grows rather rapidly, by ~ 1.4GeV~^ when 
W grows by one order in magnitude from the fixed target W = 15 GeV up to the HERA 
energy W = 150 GeV. This corresponds to the effective shrinkage rate a'^j^j ~ 0.15 GeV~^, 
only at very high energy beyond the HERA range the shrinkage will follow the true slope 
of the Regge trajectory for the rightmost gBFKL singularity a'jp = 0.07 GeV^^. The values 
of a'jp and a'^j^j- are very sensitive to the gluon propagation radius Rc and can eventually 
be used to fix this poorly known parameter. At the moment we explore major properties 
of the solution for Rc = 0.27fm. In Fig. 11 we also show the diffraction slope for the J/\E' 
production at = 100 GeV^ which nearly coincides with that for real photoproduction of 
the T. This is still another example of flavor symmetry restoration, because the scanning 
radii rs for the two reactions are very close to each other. 

The flavor symmetry properties of the diffraction cone can be seen in Fig. 12. The curves 
for -6(7* —>■ V) of all the vector mesons do converge together as a function of Q'^ + rriy, slight 
departures from exact flavor symmetry can be attributed to a difference of the pQCD scale 
factors t{V) for light and heavy vector mesons. At fixed W, the calculated dependence is 
an interplay of the changing scanning radius rs and of the decrease of the Regge component 
with the increase of Xg//. Fig. 4 shows that for the light vector mesons and < 10 GeV^ 
the soft pomeron is substantial and the high precision experimental data on the p^, (fP in 
this region of are indispensable for better understanding of the soft pomeron. In Fig. 12b 
the same results are presented as a function of the scanning radius rs as defined by Eq. (0) 
with A = Q. Crude estimates for the dependence of -8(7* V) reported in [EH] are close 
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to the present results. 

We can suggest useful empirical parameterizations for the diffraction slope. For pro- 
duction of heavy quarkonia, V = J/\E', T, the dependence of the diffraction slope at 
W = 100 GeV and in the considered range of < 500 GeV^ can be approximated by 

5(7* ^V)^Po~Pi log (91±I!±\ (64) 

with the slope /3i ^ 1.1 GeV"^ and the constant /Sq ^ 5.8GeV~^. Although (|6^) must be 
regarded only as a purely empirical crude parameterization, the logarithmic term ( |6^ is 
suggestive of a substantial role of the term (^D|) in the diffraction slope at high energy. We 
recall that the constant is subject to the choice of the t range, it is value of the slope 
/?! which is more closely related to the gBFKL djTiamics. For the light vector mesons, a 
somewhat better approximation to the results shown in Fig. 12 is 

with the same (3i = 1.1 GeV^^ as above and with (3q = 7.1 GeV~^, (32 = 1-6 for the p° 
production and j3o = 7.0 GeV~^, = 1-1 for the 0° production. 

The experimental studies of the and energy dependence of the diffraction slope are in 
the formative stage. In the heavy quarkonium sector, only photoproduction of the J/'^ has 
been studied to some extent. The experimental situation is summarized in Fig. 13, at the 



both fixed target p2| , |5T| , |61[] and HERA energy [^, ^ |5^, ^ the error bars are too big 
for the definitive conclusions on the presence and/or lack of the shrinkage of the diffraction 
cone to be drawn. On the experimental side, the determinations of the diffraction slope 
are very sensitive to the rejection of the inelastic background. Only the E401 experiment 



51| has used a technique which allowed a direct selection of the purely elastic events. The 
E401 result B{W = 15 GeV, = 0) = 5.6 ± 1.2 GeV^^ jg consistent with the NMC result 
B{W = 15 GeV, g2 = 0) = 5.0 ± 1.1 GeV^^ |6|. The recent high statistics Fermilab E687 
experiment has used the nuclear target and its determination of the diffraction slope 
for the quasielastic scattering, B{W = 20 GeV, = 0) = 7.99 ± 0.23 GeV'^, is subject to 
the model-dependent separation of the coherent and quasielastic production on nuclei. At 
HERA, the first HI data gave B{W = 90 GeV, = o) = 4.7 ± 1.9 GeV^^ [|6| and the 
first ZEUS data gave B{W = 90 GeV, Q"^ = 0) = 4.5 ± 1.4 GeV'^ |3|, updated with higher 
statistics to B{W = 90 GeV, = 0) = 4.6 ± 0.6 GeV^^ [g. In 1996 the HI collaboration 

57| found weak evidence for shrinkage of the diffraction cone: B{W ~ 60 GeV, = 0) = 



3.7 ± 0.2 ± 0.2 GeV~2 _ ^20 GeV, = 0) = 4.6 ± 0.3 ± 0.3 GeV-^ 
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For virtual production of J/\I/ the HI |£0[ reported in 1996 the first data: B{W = 
90 Gel/, < g2 >= iSGeV^) = 3.8 ± 1.2(st)t?:[j(syst) GeV^^^ Recently the ZEUS collabo- 
ration |59[ presented the value of the diffraction slope at =6 GeV^, B{W = 90GeV, < 
q2 >=6GeV^) = 4.5±0.8(st)±1.0(s2/st) GeV"2. We predict the decrease of the diffraction 
slope from Q"^ = to Q"^ = 18 GeV^ by mere ^ 0.5 GeV~^, too small an effect to be seen at 
the present experimental accuracy. 

The end-point contribution from asymmetric large size dipoles with z{l — z) <^ 1 is 
different for the production of the T and L polarized vector mesons and makes the average 
scanning radius somewhat larger in the case of the T polarization. Consequently, one 
would expect the inequality of diffraction slopes Bx > Bi for the polarization states. This 
is indeed the case, but the effect shown in Fig. 14 is negligible numerically even for the 
charmonium states, because in the nonrelativistic quarkonium the end-point effects are 
strongly suppressed. For the bottonium states the Bt — Bl is absolutely negligible. 

9 What is special about diffractive production of the 
radially excited states V{2S)? 

In the diffraction production of radially excited 25* states (\1'', T') the keyword is the node 
effect. The radial wave function of the 25* state changes the sign at r ~ Rv{lS) and 
there are cancelations of contributions to the production amplitude from large size dipoles, 
r > Ry{1S), and small dipoles, r < Rv{lS), which were noticed for the first time in 
0, ^. Manifestations of the node effect for diffractive production of light vector mesons off 



nuclei have been discussed in ||30|, ^ . The detailed analysis of the forward real and virtual 
photoproduction of light 25* states (p',0') at high energies is presented in |TD[. The major 
subject of the present discussion is new manifestations of the node for the diffraction cone. 

The cancellation pattern obviously depends on the relationship between rs and position 
of the node r„ which is close to the radius of the 15* state, r„ ~ _Rv'(15'). If rs <^ -Ry(lS'), 
then the wrong-sign contribution to the production amplitudes from dipoles with r > r„ 
is small and cancelations are weak (the undercompensation scenario of Ref. [0). If 
'^s' ^ Rvi^S), then the production amphtude can even be dominated by the wrong-sign 
contribution from r above the node position (the overcompensation scenario). For the 
heavy quarkonia the scanning radius rs is sufficiently small even at = and the under- 
compensation scenario is realized. 

At fixed target energies, the node effect is sufficiently strong suppresses the ratio R21 (t = 
0) = da{'$')/da{J/'$)\t=o by almost one order in magnitude (Fig. 15). Evidently, the smaller 
is the scanning radius the smaller is the large-r contribution and the weaker is the node 
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effect, so that the ratio da{^') / da{J /'^)\t=Q rises with as shown in Fig. 15. When the 
node effect is strong which is the case for the at = 0, then even shght variations of 
the scanning radius rs can change the strength of the node effect substantially. For this 
reason one must not be surprised that at fixed target energies the ratio dcT{^') / da{J /'^)\t=Q 
changes with quite rapidly, on a scale of smaller than the natural scale my. The 
predicted energy dependence of da{^') / da{J /^)\t=Q derives from the faster growth with 
energy of the dipole cross section for smaller dipoles which also diminishes the node effect. 
In Fig. 16 we show in more detail for the HERA energy W = 100 GeV the dependence 
of the ratio of the t integrated cross sections a{2S)/(T{lS) evaluated using the diffraction 
slope B{2S) described below. The predicted and W dependence of the node effect is 
sufficiently strong to be observed at HERA. Because for the heavier h quarks the scanning 
radius in comparison to Ry{lS) is substantially smaller than for the charmed quarks, the 
node effect in the bottonium production is much weaker, the ratio d(j{T') / da{T) \t=o is larger 
and exhibits much weaker and W dependence than for the charmonium states (Fig. 15). 

The node effect is slightly different for the T and L polarization. This is best seen in 
Fig. 6 which shows that the ratio RltC^S) for the ^^'(25) production which is different from 
Rlt{^S) both in the magnitude and Q^, W dependence. 

The new effect which we focus here on is a nontrivial impact of the node effect on the 
diffraction cone. In the conventional situation the larger are the radii of the participating 
particles the larger is the diffraction slope and for real photoproduction we have a clear 
hierarchy 

5(7 ^ p°) > 5(7 ^ 0°) > 5(7 ^ J/^) > 5(7 ^ T) (66) 

which follows the hierarchy of radii RpO > i^^o > Rj/^ > Rr- Although the mean squared 
radius of the is about twice as large than Rj/q,,, the color dipole approach uniquely 
predict 5(7 \E'') < 5(7 J/^) in a striking defiance of the hierarchy (|66D, which 
has the following origin: Let M.^ and 7W> be the moduli of contributions to the V{2S) 
production amplitude from color dipoles with size r below and above the position of the 
node in the 25" radial wave function and let 5< and 5> be the diffraction slopes for the 
corresponding contributions. Because of the hierarchy (|66D we have a strong inequality 

5> > 5< . (67) 

For production of the ^(15') state 5(15) ~ 5<. Now, the total V{2S) production amplitude 
equals M.{2S) = A^< — M.^ and for the diffraction slope we find 
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which gives an estimate 

B{2S) - B{1S) ^ -(S> - B<) ^^_^^^ < . (69) 

The weaker is the node effect the smaller is the difference of diffraction slopes 5(25')— -8(15). 
The typical color dipole sizes r which enter M< and M> differ by ~ Rv{lS) and the 
difference of slopes By — B^ can be evaluated as a variation of the diffraction slope -8(15') 
when the scanning radius rs changes by the factor ~ 2 from it's value at = 0. Then 
the parameterization gives for heavy quarkonia an estimate i?> — i?< ~(l-2)/3i ~ IGeV"^. 
Eq. (^) shows the splitting B{2S) — B{1S) is further suppressed if the node effect is weak, 
i.e., if M> < M<. 

The results for the 5(15) — -6(25) are presented in Fig. 17. For the bottonium family the 
node effect is negligibly weak, cf. Fig. 15, whereas for the charmonium family the chances 
of the experimental observation of the inequality -6(25) < -6(15) are not nil at least in real 
photoproduction and in the fixed target experiments. The difference of diffraction slopes 
-6(15) — -6(25) is larger for the L polarization conforming to a stronger node effect as seen 
in Fig. 15. As discussed above and shown in Fig. 15 the node effect diminishes with energy 
and the difference of diffraction slopes -6(15') —-6(25') drops by the factor ~ 2 from the fixed 
target to HERA energies. It vanishes at large following the demise of the node effect in 
Fig. 15, the remarks on the rapid variation of the node effect on a scale of smaller than 
rriy at fixed target energies ~ 15 GeV are fully relevant to -6(15) — -6(25) too. 

Another demonstration of the node effect leading to inequality, -6(25') < -6(15'), is 

presented in Fig. 18 in the form of the t- dependence of the differential cross section ratio 

da{^ _ V')/dt 
^^''^^^^ - da{^ - V)/dt 

for real photoproduction. The ratio -R*'/(j/-^)(t) rises with t a.tW = 15 GeV as a consequence 

of the node effect, whereas atW = 100 GeV Fig. 18 shows practically constant t- dependence 

of -R-i"/(j/-i')(t) and Rr'/r(t) corresponding to a weaker node effect at larger energy and for 

heavier vector mesons, respectively, see also Fig. 17. 

There is a solid experimental evidence for the node effect in real photoproduction of the 

In 1996 the HI collaboration reported the first observation of real photoproduction of the 

^' at HERA with the result R21 = (1(7 ^ ^')/^(7 ^ J/i') = 0.15 ± 0.05 |3|. Concerning 

the fixed target experiments, only E401 has used the hydrogen target with the result a{'~f — > 

^'')/cr(7 J/^p) = 0.20 ± 0.05 at 1^ = 17 GeV. Nuclear targets have been used in all 

other experiments. Evaluation of the cross section ratio for the hydrogen target from these 

data requires corrections for the nuclear shadowing of the J/\l/ and nuclear antishadowing 

of the production, there are also systematic uncertainties with the separation of coherent 



28 



and incoherent production. Specifically, within the same color dipole model as used in this 
paper it has been shown that the ratio R21 = cr(7 — > \E'')/(t(7 J/i') is enhanced 
in incoherent production off nuclei by the factor 1.26, 1.55 and 2.16 for the Be, Fe and 
Pb nuclei, respectively. For the relatively dilute ^Li target the enhancement factor can be 
estimated as ^ 1.1. Then, the E687 result R2iiEQ87) = 0.21 ± 0.02 for the Be target at 
W = 19GeV g corresponds to R2i{E687;N) = 0.17 ± 0.02 for the free nucleon target. 



the NMC resuh R21 = 0.20 ± Om{stat.) ± 0.07{syst.) [gT| after correction for the last 
value [ID of the branching ratio BR{J/'^ Ai^/U") = 5.97 ± 0.25% gives R21 = 0.17 ± 
0.04:{stat.)±0.04:{syst.) for the passive concrete absorber atW = 15GeV which corresponds 
to R2i{NMC; N) = 0.13 ± 0.05 for the free nucleon target. The NA14 resuh 0.22 ± 0.05 
for the Li target at W = 14GeV g corresponds to R2i{NAU;N) = 0.2 ± 0.05 for the 
free nucleon target and the SLAG result 0.22 ± 0.08 for the Be target at W ^ 6.5 GeV 
||65| corresponds to R2i{SLAC] N) = 0.18 ± 0.07 for the free nucleon target. In Fig. 19 
we compare our prediction for i?2i(A^) = cr(7 —>■ ^')/cr(7 — * J/'ip) for real photoproduction 
on protons with the HI and E401 data for the proton target and the above evaluations of 
R2i{N) from the nuclear target data. The overall agreement between theory and experiment 
is satisfactory. In view of the steady collection of the data at HERA, the higher precision 
fixed target data on the hydrogen target are highly desirable to check unambiguously the 
predicted rise of R2i{N) with energy. 

10 The summary and conclusions 

The major focus of this work has been on the forward cone for diffractive real and virtual 
photoproduction of ground (IS) and radially excited (25*) states of heavy quarkonia in the 
framework of color dipole running gBFKL approach. We presented a detailed discussion of 
the color dipole factorization for diffractive amplitudes and of the relevant pQGD factor- 
ization scales with a strong emphasis on restoration of the flavor symmetry in the variable 
+ rriy. We based our analysis on solutions of the gBFKL equations for the dipole cross 
section ||26|, |[ and for the diffraction slope for the color dipole scattering amphtude p8[ . 



Starting from the same dipole cross section which provides a good quantitative descrip- 
tion of the rise of the proton structure function at small x we found encouraging 
agreement with the experimental data on the and energy dependence of diffractive J/\E' 
production. There are many interesting predictions for the J/\E' production to be tested, 
for instance, the dependence of the effective intercept Ag//. 

A detailed analysis of the energy and dependence of the diffraction slope -8(7* V) 
for vector meson production is presented here for the first time. Of primary interest is 
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the shrinkage of the diffraction cone which follows from the finding that the gBFKL 
pomeron is a set of moving poles. We identified different sources of the dependence 
of the diffraction slope. Based on the solution of the running gBFKL equation for 
the diffraction slope, we presented detailed calculations of the and W dependence of 
5(7* — > V). The present experimental data on -8(7* for the J/\E' production are not 

yet accurate enough to rule in or rule out our predictions for the shrinkage of the diffraction 
cone. 

Diffractive production of the radially 25" mesons T') is a subject on its own. The 
key new feature of production of the 25" states is the node effect, the destructive interference 
of contributions to production amplitude from small and large color dipoles because of the 
node in the radial wave function of 25* radial excitations. The resulting strong suppres- 
sion of the photoproduction agrees with the available experimental data. An interesting 
prediction from the color dipole dynamics to be tested is a rise of the cross section ratio 
(t(7 ^')/o"(7 ^/^) by the factor two from the CENR-FNAL to HERA energies. The 
new consequence of the node effect which we discussed in this paper is a counterintuitive 
inequality of diffraction slopes 5(7 \E'') < -8(7 J/^) to be contrasted to a familiar 
rise of the diffraction slope for elastic scattering processes with the rise of the radius of 
the beam and target particles. The scanning phenomenon allows to control the node effect 
varying the scanning radius with and we present the corresponding predictions for the 

dependence of the cross section ratio (7(7 \1'')/(t(7 J/'^) and of the difference of 
diffraction slopes 5(7 < B{'~f J/^)- The predicted effects for the charmonium 

family are within the reach of modern experiments. The present analysis of diffractive pro- 
duction of heavy mesons provides a useful benchmark for future applications to light vector 
mesons. The experimental comparison of virtual and real photoproduction of vector mesons 
will shed light on the transition between the soft pomeron exchange which dominates for 
the p^,uj^,(j)^ production at small and moderate to the gBKFL pomeron exchange at 
higher and/or heavy vector mesons. 
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Figure captions: 

Fig. 1 - The perturbative QCD diagrams for vector meson production. 

Fig. 2 - Tlie dipole size dependence of tlie diffraction slope for the perturbative gBFKL, 
Bpt{r), and soft, Bnptij), pomerons. 

Fig. 3 - The color dipole model predictions for the total cross section atotiVN) for the 
interaction of the heavy vector mesons J/\E',\E'', T and T' with the nucleon target 
as a function of c.m.s. energy W . The dashed shows represent the soft pomeron 
contribution. The right box shows the color dipole model predictions for the total 
cross section atotiVN) vs. the radius Ry of vector mesons p°, p', 0°, 0', J/^ ,T and 

Fig. 4 - Decomposition of production amplitude for longitudinally polarized vector mesons 
into soft (dashed curves) and perturbative+soft (solid curves) contribution as a func- 
tion of my + Q^. The non-monotonic dependence is due to the increase of Xe// at 
fixed W. 

Fig. 5 - The color dipole model predictions for the W dependence of the real photoproduction 
cross section o"(7* — > V) for the J/\E' production vs. the low-energy EMC E516 
13, E401 E687 m and high-energy ZEUS ||, ||, |5|] and HI H, |3, m data. 



Fig. 6 - The color dipole model predictions for the and W dependence of the ratio of 
the longitudinal and transverse differential cross sections in the form of the quantity 
Rlt = '^ daT^jl'^v] ' '^here mv is the mass of the vector meson. 

Fig. 7 - The color dipole model predictions for the polarization-unseparated forward differ- 
ential cross section (top boxes) (icr(7* — > V)/dt\t=o = daxij* V) / dt\t=Q + daLij* 
V)/dt\t=Q for the J/\E' and T production as a function of the c.m.s. energy W at differ- 
ent values of . The bottom boxes are predictions for the polarization-unseparated 
t-integrated cross section a{j* —^V)= aripf* ^ V^) + (7^(7* — ^ V). 

Fig. 8 - The dependence of the effective intercept /S.effiQ'^) for the forward production 
of p°, 0°, J/^ and T. 

Fig. 9 - The color dipole model predictions for the dependence on the scaling variable + 
of the polarization-unseparated da{j* — » V)/dt\t=Q = dax^Ci* — * V)/dt\t=o + 
daL{l* V)/dt\t=o at the HERA energy W = 100 Gel/. 
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Fig. 10 - The color dipole model predictions for the dependence of the observed cross 
section (1(7* V) = (Jt{^* — > V) +e(JL(7* of exclusive J/\E' production vs. the 

low-energy (EMC [||]) and high-energy (ZEUS ||, HI |8|, 0) data. 



Fig. 11 - The color dipole model predictions for the cm. s. energy dependence of the diffraction 
slope for real photoproduction of the J/\E' and T and for the J/^ electroproduction 
at g2 = lOOGeV^. 

Fig. 12 - The color dipole model predictions for the diffraction slope in production of different 
vector mesons as a function of the scaling variable my-|-Q^ (left box) and the scanning 
radius (right box) at fixed c.m.s. energy W = 100 GeV. The scales of on the top 
of the right box show the values of which correspond to the scanning radii shown 
on the bottom axis. 

Fig. 13 - Comparison of the color dipole model prediction for c.m.s. energy W dependence 
of the diffraction slope for photoproduction of the J/\l/ with the E401 NMC ^1 



E687 |32|, HI pi 13 and ZEUS |m EB data. 



Fig. 14 - The color dipole model predictions for the difference of diffraction slopes Bt — Bl 
as a function of for the J /'^ and ^' production. 

Fig. 15 - The color dipole model predictions for the and W dependence of the ratios 
da{Y ^ '^'{2S))/da{n* J/^) and c/a(7* ^ V {2S))/ dain* ^ T^) for the T, L and 
polarization-unseparated forward differential cross sections. 

Fig. 16 - The color dipole model predictions for the dependence of the ratio of the t- 
integrated polarization-unseparated cross sections cr(7* \E''(2S'))/cr(7* J/^) and 
(t(7* T'(2^))/ct(7* T) at c.m.s. energy W = 100 GeV. 

Fig. 17 - The color dipole model predictions for the difference of diffraction slopes -8(15*) — 
B{2S) vs. at c.m.s. energy W = 15 GeV (dashed lines) and W = 100 GeV (sohd 
lines) for T and L polarization. 

Fig. 18 -The color dipole model predictions for t- dependence of the ratio Ry'/vit) = '^l^l^y^^j'ft ■. 
for the $'/(J/^) and T'/T real photoproduction. 

Fig. 19 - Comparison of the color dipole model prediction for c.m.s. energy W dependence 
of the ratio ^(7 ^')/^(7 ^ J/'^) at = q with the E401 H, NMC ll, E687 
[H, NAM H, SLAG H and HI H data. 
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